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SUMMARY 

A three-bay, space, cantilever truss is probabilistically evaluated using the 
computer code NESSUS (Numerical Evaluation of Stochastic Structures Under Stress) 
to identify and quantify the uncertainties and respective sensitivities associated 
with corresponding uncertainties in the primitive variables (structural, material, 
and loads parameters) that defines the truss. The distribution of each of these 
primitive variables is described in terms of one of several available distribu- 
tions, such as the Weibull, exponential, normal, log-normal, etc. The CDF’s for 
the response functions considered and sensitivities associated with the primitive 
variables for given response are investigated. These sensitivities help in deter- 
mining the dominating primitive variables for that response. 

INTRODUCTION 

It is a common practice to evaluate the structural integrity of trusses with 
the aid of deterministic analysis techniques and appropriate safety factors. Tra- 
ditionally, these factors are an outcome of many years of analytical as well as 
experimental experience in the area of structural mechanics design. The load fac- 
tors are used to take into account many different operating conditions and also to 
ensure that the maximum operating load has been considered. The safety factors are 
used to account for unknown effects in analysis and fabrication. 

As an alternative to the deterministic approach, probabilistic structural 
analysis methods (PSAM) have been developed at NASA Lewis Research Center (ref. 1) 
which consider various uncertainties in a more structured manner. PSAM takes into 
account the scatter in the resistance (strength) and parameters such as geometric 
configuration, loadings, structural properties, etc., herein referred to as primi- 
tive variables. Furthermore, the primitive variables are not considered as either 
single values or upper and lower bound values. Instead, the actual probability 
distribution of the primitive variables is represented. Thus, the probabilistic 
methodologies have taken a prominent role in designing complex structural compo- 
nents which have experienced failures using a deterministic analysis approach 
(ref. 2). 

The NESSUS Computer Code developed at NASA Lewis (refs. 3 and 4) contains PSAM 
techniques and provides a choice of algorithms for the solution of static, dynamic, 
buckling, and nonlinear analysis problems. It also includes a number of innovative 
analyses to evaluate the sensitivity of response variables to small variations in 
one or more user-defined primitive variables. Recently, NESSUS has been used for 
the analysis of Space Shuttle Main Engine (SSME) components, for example, to obtain 
a probabilistic assessment of a mistuned bladed disk assembly (ref. 5) and to eval- 
uate the reliability and risk of a turbine blade under complex service environments 
(ref. 6). 

The objective of this paper is to demonstrate the use of the NESSUS Computer 
Code to probabilistically evaluate a three-dimensional, three-bay, cantilever truss 



typical for space station structures with the aim to identify and quantify the sen- 
sitivities associated with uncertainties in primitive variables. For a space truss 
the primitive variables such as spatial truss geometry, stiffness parameters, 
strength parameters, and applied loads or moments continuously vary due to changes 
in the environment. The distribution of each of these primitive variables is 
described in terms of any one of several available distribution functions, such as 
the Weibull, exponential, normal, log-normal, etc. These distributions have sig- 
nificant impact bn the scatter of the response variables such as nodal displace- 
ments, member forces, vibration frequencies, etc. The specific PSAM technique used 
in the NESSUS computer code evaluates the scatter of the response variables from 
the scatter in the primitive variables. 


PROGRAM CAPABILITY AND DESCRIPTION 

The NESSUS code consists of three major modules namely; NESSUS/PRE, NESSUS/ 
FEM, and NESSUS/FPI. It is important to note that each of these modules can be 
used independently. NESSUS/FEM and NESSUS/FPI are combined into a single computer 
code called NESSUS/PFEM which performs the entire probabilistic finite element 
analysis including perturbations of the primitive variables. In general, the prim- 
itive variables are specified with their mean values (n), standard deviation fa) 
and the type of distribution. 


The NESSUS/PRE module is a pre-processor used to obtain the description of a 
partially correlated Gaussian field in terms of a set of uncorrelated random vec- 
tors. The NESSUS/FEM module is a finite element analysis code with the capability 
to generate perturbed solutions about a deterministic state. It contains an effi- 
cient perturbation technique such that the perturbation of each variable is done 
rapidly. Each perturbation corresponds to a prescribed deviation from the determi- 
nistic model. The NESSUS/FPI module contains several advanced reliability methods 
including Monte-Carlo simulation. The fast probability integration (FPI) (ref 6) 
techniques are one or several orders of magnitude more efficient than the Monte 
r meth ? ds ; The module extracts the database of perturbed solutions 

from NESSUS/FEM to calculate the probability distribution functions of the response 
variables. 


PROBABILISTIC FINITE ELEMENT ANALYSIS 

Dias et al. (ref. 3) discuss the developments in probabilistic finite element 
analysis which are based upon perturbation methods for computing the sensitivity 
of the response of the random variables present. In general, the finite element 
equation for motion is written as: 


[M]{u} + [C]{6} + [K] {u} = F(t) ( 1 ) 

where [M], [C] and [K] denote the mass, damping, and stiffness matrices, respec- 
tively. It is important to note that these matrices are calculated probabi 1 isti- 
cally in the NESSUS code. Furthermore, {u} , {u} , and {u} are the acceleration, 
velocity and displacement vectors at each node, respectively. The forcing function 
vector, {F ( t ) } , is time dependent at each node. 

In this paper, the static case is considered by setting the mass and damping 
matrices to zero and considering the forcing function being independent of time in 
equation (1) such that 


[K]{u} = {F} 


( 2 ) 



Furthermore, the eigenvalue analysis is also carried by setting only the damping 
matrix to zero and using the following equation; 

[K] - w 2 [M] {u} = 0 (3) 

where w denote eigenvalues and {u} are corresponding eigenvectors. 


FINITE ELEMENT MODEL 

A three-dimensional, three-bay cantilever truss is simulated using a linear 
isoparametric beam element based on the Timoshenko beam equations. The element is 
idealized as a two-noded line segment in three-dimensional space. The cantilever 
truss is assumed to be made from 44 hollow circular pipe type beam elements (see 
fig. 1). The pipes are made up of wrought Aluminum alloy (616-w) with modulus of 
elasticity (E) equal to 10 mpsi. The outer and inner radii (r 0 and ri) of the 
tube, are 0.5 and 0.4375 in., respectively. All 6 degrees of freedom are 
restrained at the fixed end (left side) nodes. The truss is analyzed twice, once 
using beam elements and then using pseudo-truss elements. The beam element is con- 
verted into a pseudo-truss element by suppressing the effective shear areas in the 
principal planes (A xx and Ayy), the two principal moments of inertias for the 
tube cross-section, (I xx and Iyy) and torsional constant, J. In the case of 
truss elements, 3 rotational degrees of freedom at each node and 3 translational 
degrees of freedom at support nodes are restrained. 

Each bay of the truss is 5 ft wide, 8 ft long, and 6 ft high (fig. 1). The 
overall length of the truss is 24 ft. Six vertical and two longitudinal loads are 
applied. Twisting moments are applied at the truss-end top nodes with beam ele- 
ments. However, an equivalent couple is applied at the truss-end nodes for truss 
elements. The directions of the forces and moments are shown in figure 1 and the 
mean values are given in table I. 


PROBABILISTIC MODEL 

The following primitive variables are considered in perturbation analysis: 

(1) Nodal Coordinates (X,Y,Z) 

(2) Modulus of elasticity (E) 

(3) Outer radius of the tube (r 0 ) 

(4) Inner radius of the tube (rj) 

(5) Vertical loads (V) 

(6) Longitudinal loads (H) 

(7) Truss-end moments (M) 

(8) Truss-end coupling forces (P) 

It is possible that the above design variables will vary continuously and 
simultaneously due to extreme changes in the environment when such trusses are used 
in upper earth orbit for space station type structures. The normal distribution is 
used to represent the scatter in E, r 0 , rj, and X,Y,Z coordinates. The applied 
loads, moments and coupling forces are selected to represent anticipated loading 
conditions for a typical space truss. These are represented by log-normal distri- 
butions. Initially, the deterministic finite element analysis takes into consider- 
ation the mean value of these primitive variables. In the probabilistic analysis 
each variable is perturbed equidistant from the mean value. However, each variable 
is perturbed independently and by a different amount. Usually, the perturbed value 
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of the design variable is obtained by certain factor of the standard deviation at 
either side of the mean value. Finally, the NESSUS/FPI module extracts response 
variable values (one deterministic and two times the number of primitive variables) 
to calculate a probability distribution function of the response variable consid- 
ered. The mean, distribution type and percentage variation for different primitive 
variables are given in table I. 


DISCUSSION OF RESULTS 

The three-bay cantilever truss is probabilistically analyzed and the cumula- 
tive probability distributions for the truss end displacements, member forces and 
vibration frequencies are plotted. The sensitivities of the primitive variables on 
the scatter in the truss structural responses (truss free end displacements, member 
axial forces and vibration frequencies) are quantified in table II. Figures 2 to 4 
depict the probabilistic displacement of the truss free end nodes (top and bottom) 
in X, Y, and Z directions, respectively, using the truss element. The large dif- 
ferences in CDF's between top and bottom nodes as seen in figures 3 and 4 are, 
respectively, due to longitudinal loads at top nodes and vertical loads at bottom 
nodes. However, the results obtained using beam elements show a lower magnitude 
of the probabilistic displacement in Z direction as seen from figure 5. This is 
because of the fact that the beam elements increase the over all stiffness of the 
truss. It is noted from the table II that the perturbations in Y coordinates have 
a major impact on the scatter in the displacements. In addition, the perturbations 
in Z coordinates have an impact on truss end displacements in the Z direction when 
beam elements are used. The CDF of the axial forces in the Top Front Longeron 
(TFL) , Bottom Batten (BB), Front Vertical (FV), and Rear Diagonal (RD) of the first 
bay from the support are shown in figures 6 to 9. The changes in axial forces, 
from compression to tension, in Bottom Batten (fig. 7) and Front Vertical (fig. 8) 
are due to reversal in loading or due to loading from either side of the mean 
value. The CDF's of frequencies of modes 1 and 2 using truss elements are plotted 
in figures 10 and 11, whereas, figures 12 and 13 are those for the truss with beam 
elements. It is important to note from table II that the cross sectional area 
(primitive variables r 0 and r^) has a significant impact on the probabilistic 
distribution of the vibration frequencies. Furthermore, the truss with beam ele- 
ments has, higher frequencies due to inclusion of bending stiffness. Finally, the 
deterministic value of the response variable may be estimated from the 50 percent 
cumulative probability level. 


SUMMARY 

The application of the PSAM code to model and probabilistically evaluate a 
cantilever truss typical for space station is demonstrated using beam and truss 
elements. The probabilistic structural responses are predicted and plotted. The 
sensitivities associated with uncertainties in the primitive variables are quanti- 
fied. The results indicate substantial scatter in frequency and some member axial 
force distributions. However, there is relatively less scatter in displacement 
distributions . 
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TABLE I. - PRIMITIVE VARIABLES AND UNCERTAINTIES FOR 
PROBABILISTIC STRUCTURAL ANALYSIS OF A SPACE TRUSS 
[Random Input Data.] 


Primitive variables 

Oi stri buti on 

Mean value 

Scatter, 



type 


ipercentage 

Geometry 

Width 

Normal 

60 in. 

6.0 


Length 

Normal 

96 in. 

6.3 




192 in. 

7.5 




288 in. 

6.3 


Hei ght 

Normal 

72 in. 

7.5 

Loads 

Vertical 

Log-normal 

200 lb 

6.3 


Longi tudi nal 

Log-normal 

200 lb 

2.5 


Couple 

Log-normal 

0.7 lb 

6.3 


End moment® 

Log-normal 

50 lb/in. 

6.3 

Material 

property 

Modul us 

Normal 

10 Mpsi 

7.5 

Tube radii 

Outer radius 

Normal 

0.5 in. 

7.5 


Inner radius 

Normal 

0.44 in. 

7.5 


a For beam elements only. 
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TABLE II. - SENSITIVITIES OF PRIMITIVE VARIABLES UNCERTAINTIES OF TRUSS STRUCTURAL RESPONSE 


Response type 

Sensitivity factors 

Geometry 

Loads 

Material 

property 

Tube radii 

Width 

Length 

Hei ght 

Verti- 

cal 

Longi tu- 
di nal 

End 

moment 

Coup! e 

Modul us 

Outer 

Inner 

Di spl acement: 
X-di recti on 
Y-di recti on 
Z-direction 
Z— d i recti on a 

0.35 

.18 

.12 

(b) 

0.66 

.69 

.73 

.76 

(b) 

(b) 

(b) 

0.34 

0.22 

.36 

.28 

.34 

(b) 

0.16 

(b) 

(b) 

(t 

>) 

0 

>) 

(t 

>) 

D 

(b) 

Axial force: 

Top front longeron 
Bottom batten 
Front vertical 
Rear diagonal 

0.24 

.39 

(b) 

.55 

0.44 

.28 

(b) 

.39 

(b) 

0.14 

(b) 

(b) 

II 

HI 

(1 

>) 

(1 

>) 

(l 

>) 

(b) 

(b) 

0.12 

(b> 

i 

(b) 

0.44 

.63 

(b) 

Frequency: 
Mode-1 
Mode-1 a 
Mode-2 
Mode-2 a 

1 

0.47 

.30 

.39 

.27 

(b) 

0.11 

(b) 

.18 

(b) 

<t 

>) 

d 

>) 

O 

>) 

* 

(t 

>) 

0.25 

.30 

.32 

.33 

0.19 

.23 

.24 

.25 


a For beam element only. 

“Sensitivity factors less than 10-percent. 



FIGURE 1. - SOLAR ARRAY PANELS MAST - TYPICAL TRUSS. 


FIGURE 2. - PROBABILISTIC DISPLACEMENT. 
TRUSS FREE END: X - DIRECTION' (WIDTH). 


6 













































CUMULATIVE PROBABILITY CUMULATIVE PROBABILITY 



> 

_J 

CD 

< 

CQ 

0 

1 

a 


LU 

> 

R 

< 

D 

2 

=> 



FIGURE 3. - PROBABILISTIC DISPLACEMENT. 
TRUSS FREE END: Y - DIRECTION (LENGTH), 


FIGURE A. - PROBABILISTIC DISPLACEMENT. 
TRUSS FREE END: Z - DIRECTION (HEIGHT). 



DISPLACEMENT (in.) 

FIGURE 5. - PROBABILISTIC DISPLACEMENT. 
TRUSS END: Z - DIRECTION (HEIGHT). 



FIGURE 6. - PROBABILISTIC MEMBER FORCE. 
FIRST BAY TOP FRONT LONGERON. 
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CUMULATIVE PROBABILITY CUMULATIVE PROBABILITY 



AXIAL FORCE (lbs) 

FIGURE 7. - PROBABILISTIC MEMBER FORCE. 
FIRST BAY RIGHT BOTTOM BATTEN. 


- 505.50 - 153.30 198.90 551.10 903.30 

AXIAL FORCE (lbs) 

FIGURE 8. - PROBABILISTIC MEMBER FORCE. 
FIRST BAY RIGHT FRONT VERTICAL. 



AXIAL FORCE (lbs) 

FIGURE 9. - PROBABILISTIC MEMBER FORCE. 
FIRST BAY REAR DIAGONAL. 



FIGURE 10. - PROBABILISTIC FREQUENCY MODE 


1255.50 


18.00 





FREQUENCY (cps) 

FIGURE 11, - PROBABILISTIC FREQUENCY MODE 2. 




FREQUENCY (cps) 

FIGURE 13. - PROBABILISTIC FREQUENCY NODE 2, 
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